Stationary Configurations Imply Shift Symmetry: No Bondi Accretion for 

Quintessence / k-Essence. 
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In this paper we show that, for general scalar fields, stationary configurations are possible for 
shift symmetric theories only. This symmetry with respect to constant translations in field space 
should either be manifest in the original field variables or reveal itself after an appropriate field 
redefinition. In particular this result implies that neither k-Essence nor Quintessence can have exact 
steady state / Bondi accretion onto Black Holes. We also discuss the role of field redefinitions in k- 
Essence theories. Here we study the transformation properties of observables and other variables in 
k-Essence and emphasize which of them are covariant under field redefinitions. Finally we find that 
stationary field configurations are necessarily linear in Killing time, provided that shift symmetry 
is realized in terms of these field variables. 



INTRODUCTION 



The discovery of the late time acceleration of the uni- 
verse using Supernova la [3, [2] confirmed by other ob- 
servations (see Ref. 0] and references therein), opened a 
window of opportunity for the existence of novel cosmo- 
logical scalar fields not only during the early inflationary 
stage but also in the current universe. Indeed the scalar 
fields are the most natural candidates for realization of 
inflation and for the dynamical explanation of Dark En- 
ergy (DE) which is responsible for the late time acceler- 
ation. Arguably, the main difficulty in the modeling and 
understanding of the possible dynamics of Dark Energy, 
arises because of the fine tuning issues. Inparticular, 
there is the so-called coincidence problem 0, 0| : why the 
energy density in DE is only now comparable with the 
energy density in the dust-like Dark Matter? This coin- 
cidence would be especially remarkable, if one assumes 
that both these Dark constituents are independent of 
each other and evolve very differently in time. Partially 
because of the fine tuning problems it is not surprising 
that the candidates for DE often have not only rather 
exotic names: Quintessence/ Cosmon [El. IfjL la. IsL [Io| . 
k-Essence [iH fl2l Il3l ] . Phantom Ghost Condensate 
[l5| . Quintom |la | etc but also correspondingly very un- 
usual properties. In particular, these scalar fields can 
possess: extremely small effective mass (Quintessence, 
Quintom), sound speed which can be much smaller and 
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even larger than the speed of light (k-Essence, Ghost 
Condensate), negative kinetic energies (Phantom, Quin- 
tom), Lorentz symmetry breaking and gravity modifica- 
tions even around the Minkowski space-time background 
(Ghost Condensate). The most successful paradigm to 
solve the coincidence problem is currently the k-Essence, 
where the highly nonlinear dynamics triggers the equa- 
tion of state of DE from radiation-like to quasi de Sitter 
around the transition to the matter domination stage. In 
the late matter domination epoch the k-Essence has the 
speed of sound which is much smaller than one. How- 
ever, it was showed [TtI . [l8| that to explain the coinci- 
dence problem k-Essence models must necessarily have, 
at least, a short phase where the fluctuations in the k- 
Essence travel at superluminal speeds. For our paper it 
is important that the nonlinear dynamics responsible for 
the attractor behavior addressing the coincidence prob- 
lem requires an explicit dependence of the Lagrangian 
on the field strength [l8j]. This field dependence cannot 
be eliminated by any field redefinitions. Thus, success- 
ful k-Essence models as well as Quintessence / Cosmon 
models cannot be shift symmetric. 

On the other hand it is known that the current uni- 
verse is highly inhomogeneous on small scales and in 
particular that there are plenty of Black Holes (BHs) 
of different mass and origin. Thus an interesting and 
natural question arises, how do Black Holes surrounded 
by cosmological scalar fields evolve? In addition, from 
the theoretical viewpoint it is interesting to consider 
BHs "dressed" with different field backgrounds. This 
could have a valuable impact on our understanding of 
the physics of horizons (see e.g . figL IgcL l2ll. [2 20 ■ Ow- 



ing to the no-hair theorems [23, [2J, [25|, [26|, [27] we know 



that BHs cannot support static configurations of scalar 
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fields 1 . Therefore, any scalar hair will be continuously 
swallowed by the BH. In particular one could analyze 
the growth (and may be even formation) of Black Holes 
due to the accretion (collapse) of DE. Then one can try 
to use powerful and rather universal laws of Black Hole 
thermodynamics combined with astrophysical observa- 
tions to restrict the allowed properties of DE candidates 
and rule out some of them as contradicting to either BH 
thermodynamics or astrophysical data. Recent studies 
along these lines were, for example, done e.g. in Ref. 

0, E3, ED, m El, M, GEO, \M, \M, H EE El. 

Finally, for k-Essence, a typically very small sound 
speed during the late matter domination era allows for 
rather significant large-scale inhomogeneities around BHs 
and other massive objects. This long-range clumping 
would be one of the characteristic, potentially observ- 
able consequences of k-Essence. Moreover, due to this 
ability to realize small sound speeds along with the dust- 
like equation of state, the k-Essence fields can be used to 
model Dark Matter [13, El, EH . In this setup, the pres- 
ence of supermassive BHs at the center of galaxies makes 
understanding the accretion process even more necessary. 

On the other hand the presence of backgrounds with 
the superluminal sound speed mentioned above opens an 
exciting possibility to look beyond the BH horizon [20, 
ED] 2 . Note that the current bounds [H [H, H, [H, EH 
on DE sound speed are not restrictive at all. 

The classical and most simple setup for accretion prob- 
lems is a steady state or Bondi accretion Remark- 
ably, a lot of astrophysical phenomena can be described 
by a steady state accretion. For a review see e.g. [4gj |. 
For scalar fields, the Bondi accretion was recently studied 
in [IS ED, E3, El El SI- It is fair to say that almost all 
known analytical solutions 3 for accreting scalars either 
belong to the Bondi case or represent the dust-like free 
fall. The dust-like time dependent accretion of a mas- 
sive canonical scalar field was considered in [30], while 
dust-like solutions for the Ghost Condensate scalars were 
found in (H, . It seems that scalar fields with canon- 
ical kinetic terms would not leave any important im pac t 
on the astrophysical BHs in the current universe [301 ] - 
Nevertheless, accreting scalar fields could play an impor- 
tant role for the formation of primordial BHs (see e.g. 
0). 

In this paper we investigate stationary configurations 



1 BHs can not support scalar hair at least for theories that re- 
spect some of the standard energy conditions. Having in mind 
the exotic properties of DE models mentioned above, it would 
be interesting to find examples of stable scalar hair in theories 
violating the usual energy conditions. For a model of hairy scalar 
BHs with ghost like quantum instabilities see Ref. [28l |. 

2 Despite of the presence of the superluminal propagation the 
accretion backgrounds constructed in these works are free of 
any causal pathologies 1401 . However, it is interesting to study 
whether, similar to Ref. [4ll |. two boosted BH could create causal 
paradoxes in this setup. 

3 See however Ref. 



for general k-Essence scalar field theories. We show that 
the necessary condition for the existence of exact station- 
ary configurations is the symmetry of the theory with 
respect to constant shifts in the field space: 4> ~~ * 4> + c - 
This symmetry has to be realized either in terms of the 
original field strength or after a field redefinition. On 
the way, we also analyze properties of general k-Essence 
scalar field theories covariant with respect to field re- 
definitions. The proof is valid for general theories with 
nonlinear kinetic terms in both the test-field approxima- 
tion and the self-consistent case where the background 
metric is governed by the field 4> itself. It is interesting 
to note that shift symmetric scalar field theories are ex- 
actly equivalent to perfect fluid hydrodynamics provided 
that only such field configurations which have timelike 
derivatives are considered. In particular this result im- 
plies that the most interesting scalar field models of Dark 
Energy cannot realize a steady state / Bondi accretion. 
Thus, in general, the solution to the problem of accretion 
of these fields onto Black Holes requires a knowledge of 
their initial configuration. In this paper we are discussing 
stationary configurations, which are exact. Of course for 
the real world the stationarity should be considered as 
an approximation. It may well happen that the solu- 
tions would only asymptotically approach the stationary 
regime. For some canonical scalar fields this behavior 
was demonstrated in [30| . 



II. DERIVATION OF THE STATIONARY 
CONFIGURATIONS 

Let us consider a general scalar field theory with the 
action 

S = j d 4 aV^P (0, X) , where X = ^V^V^, 

g^ v is the gravitational metric and as usual g = detg^. 
Throughout the paper is the covariant derivative as- 
sociated with the gravitational metric g^ v . We assume 
that the Lagrangian P (0, X) is a general function sat- 
isfying the following conditions: P,x > (Null Energy 
Condition) and 2XP t xx/P,x > — 1 (Hyperbolicity con- 
dition) 4 . The first condition guaranties that the per- 
turbations carry positive kinetic energy while the second 
one implies the stability with respect to high frequency 
perturbations and is necessary for the Cau chy problem 
to be well posed (see e.g. Refs. [13, 53, H El, [Hi, E3). 
These conditions restrict the variety of the allowed La- 
grangians along with the corresponding solutions and are 
unavoidable for any physically meaningful model 5 . The 



4 In this paper we use the notation (■■■), x = d(-) /9X and the 

signature (H ). 

6 For a different opinion see 
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energy-momentum tensor of the theory is 
Tftu = FxV p 0V„0 - g^P. 



(1) 



It is well known (see e.g. [55(1), that for the timelike 
derivatives X > 0, the models under consideration can be 
described in a hydrodynamical language by introducing 
an effective four velocity 6 



V, 



along with the pressure 

P = p{4>,x), 

the energy density 

e{<j),X)=2XP tX 
and the sound speed 7 

-P. xx 



P. 



<>,X)=[1 + 2X 



dp 
Ife 



(2) 

(3) 
(4) 

(5) 



In these variables the energy momentum tensor has the 
form corresponding to the one of a perfect fluid 

Tpv = (s +p) u^Uv - pg^v 

It is convenient to use the hydrodynamical notation for 
these functions of 4> and X also for X < when they do 
not have their usual physical meaning of velocity etc. 



A. Field redefinitions and conditions for 
stationarity 

If the field <f> does not have any direct interactions 
except with gravity, then obviously a field redefinition 
<j) = <j){4>) cannot affect any observables besides the field 
itself. This is a particular case of a stronger statement 
(see e.g. Ref. [13]). Obviously the solutions (p(x) 
and 4> ( x ) result through Einstein equations in the same 
gravitational metric g^ (x) and describe in that sense 
the same physical process. Thus it is interesting to in- 
vestigate the properties of k-Essence under field redef- 
initions. Under field redefinitions 6 = 6(<b) we have 



6 Note that even for X > the effective four velocity introduced in 
J2J is not necessarily future directed. However, the analogy with 
the perfect fluid can be made exact by multiplying this expression 
J2J with ±1 so that u° > 0. Furthermore, it is convenient to use 
the analytic definition of the square root so that every time when 
4> changes its sign the square root will change the sign as well 
preserving the future direction of « M . 

7 This formula for the sound speed was introduced for the cosmo- 
logical perturbations in [5f|. One can show |iol | that the same 
expression is valid in the general case of backgrounds with time- 
like field derivatives: X > 0. 



V M = (d(f>/d<f\ V M whereas the expressions for the 

energy momentum tensor and all hydrodynamical 
quantities e, p, c s and m m remain unchanged or covariant 
8 . Here we should distinguish between covariance and in- 
variance. Covariance means that the way how the quan- 
tities / equations are constructed from other objects re- 
mains unchanged whereas invariance implies exactly the 
same functional dependence on these objects. For exam- 
ple, the formula ([4]) defining the energy density e through 
Lagrangian P, X and the derivative P t x looks the same 
after a field redefinition (covariant), however the depen- 
dence of the Lagrangian on the field does change (not 
invariant). It is obvious that e.g. the value of physical 
energy density at every point should not change under 
field redefinitions, but here these quantities reveal in ad- 
dition such covariance with respect to the field redefini- 
tions as it is the case for e.g. Euler-Lagrange equations. 
However, this covariance is not guarantied for all inter- 
esting objects. It is worthwhile mentioning that, e.g. the 
metric [Hj] 



P 



Gnu [4>o\ = I 



P 



,xx 



p 



X 



describing the propagation of small perturbations tt 
around a given background <fio (x) transforms conformally 
under field redefinitions 6 = 



G^ v [0 O ] = (jj^ G f*AM- 
Thus, as expected, the causal structure does not change 



under field redefinitions. The conformal factor 



compensates for the redefinition of perturbations n 
(dcf>/d^j Jr. 

Let us further consider a stationary space-time with 
metric g M „ and a timelike Killing vector t a . Thus 
£tg^u = 0, where £t is the Lie derivative. The con- 
figuration is stationary, if per definition 



£tPfj,f 



Using Leibniz rule we have 



0. 



£tT^ = {£ t Px) V p 0V„0 - g^£ t P + 

+Px [{£tV^) Vv<p + {£ t V„<f>) V M 0] = (6) 

By multiplying this expression with g^ v we obtain 

- £ t T£ - £ t {2XP X - 4P) = £ t (e- 3p) . (7) 

Suppose the configuration <j> (x M ) is such that V M is a 
null vector: X = 0. In that case we can multiply the right 



Note that the four velocity J2J is invariant up to the sign only. 
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hand side of the Eq. © with g^ v to obtain £ t P = 0. 
Further we have = £ t P = P^d t <f). As we are looking 
for stationary but not static solutions we have P^ — 0. 
Thus the Lagrangian should be symmetric with respect 
to field shifts <f> —> <j) + c, where c is an arbitrary constant. 
For X y^z it is convenient to introduce the projector 



rfiv — 9fJ-V 

with the properties 



2X 



(8) 



V^V<t> = 0, V^V Xv = V; and V» = 3. (9) 

Moreover, this projector is both invariant and covariant 
under field reparametrizations: V^ v [4>] = [(/)]. By act- 
ing with the projector on the left hand side of Eq. 
© we have = V fn, £ t T^ = -3£ t P. Therefore, if the 
configuration is stationary then in particular 



£ t P = 0, 



(10) 



which for the hydrodynamical case reduces to the con- 
stancy of pressure p. Combining this with (J7j) we obtain 
the time independence of the energy density e or 



£t (XPx) = 0. 



(11) 



Further we can act on the left hand side of Eq. |(6|) with 
V av : so that = V av £ t T ll , v = P x V av (£ t V v <j>) V M £ 
Thus, stationarity implies 

T a(3 {£ t V^ 



= 0. 



Using the properties of the projector ©, Leibniz rule 
and that t a is a Killing vector one obtains 







_ T> a P 



= -v, 



-pa/3 _ 



2X 



■£ t X. 



The last expression in turn can be written in the following 
form 



-£ f W c 



~2X 



£tX 



Therefore, stationarity implies 
£t ' 



= 0, 



(12) 



or in the hydrodynamical notation £ t w M = 9 . Thus 
we have proved that for any stationary configuration the 
following conditions 



£ t u li = 0, £t£ — and £tp = 0, 



(13) 



9 The vector m m is formally imaginary for X < 0. However, with- 
out any change of the results one could redefine tt M in this case: 
u* 1 = V^0/v^2X . 



should be satisfied. Note that these conditions are co- 
variant under field redefinitions and, for the hydrody- 
namical case (X > 0), are intuitively clear requirements. 
Sometimes (see e.g. [31J) one claims that the stationarity 
implies a stronger requirement: 



£Na<P = 0, 



(14) 



instead of the condition (jT2]) . However, the equation 
above is not covariant under the field redefinitions and 
does not follow from the stationarity of the energy mo- 
mentum tensor. 

Now let us find what type of theories P (</>, X) and field 
configurations <j>(x^) can, in principle, satisfy conditions 
(fl"3"ll . It is convenient to chose a coordinate system (t, x l ) 
such that the time coordinate corresponds to the integral 
curves of t a . In that case the Lie derivative reduces to 
the partial derivative £t = dt- 



B. Which filed configurations can have constant 
effective four velocity ? 

Now let us find the configurations <p (a; M ) satisfying the 
condition on the effective four velocity lfl2|) . For the time 
component of the four velocity we have 



dt 



1 



= 0, 



'2X J V2X 
where <j> — d t <j>, while for the spatial components 

di6 \ did) „ ,„ / 1 



(15) 



d, 



= 0. 



(16) 



Obviously these equations have a trivial static solution 
4> = dj^x' 1 ). To find a nontrivial solution we combine 
these two equations to obtain following system of equa- 
tions 



= 0, 



which is equivalent to 



(17) 



This is a system of partial differential equations of the 
second order. Integrating Eq. I|17|) we obtain the follow- 
ing linear homogeneous system 



did) = Vi (x 1 



(18) 



where Vi (a^) are unknown time independent functions. 
This is the first order system of three partial differen- 
tial equations for only one function <f>. Let us find the 
consistency conditions under which the system can have 
solutions. Differentiating i— equation with respect to x 1 



5 



and using the time differentiation of the j— equation we 
obtain 



djdi<j> = djViti 



Vidj4> = djVtf 



ViV j( 



Now we can compare this result with the result of the 
same procedure performed for the j— equation. We ob- 
tain 



o. 



For a simply connected manifold, the last equation im- 
plies the existence of a function (potential) \& (a; 1 ) such 
that Vi = di^. Otherwise there are no solutions for (jTSj) - 
For the i— equation we can assume that all x k with 
k ^ i are frozen parameters and for the characteristics 
(for the method of characteristics see e.g. excellent book 
581) we obtain 



dx i 
~d7 



1. 



The first integral X of this system is given by the constant 
of integration for the equation 



dt 

dx l 



= -dcS (x j 



By integrating which we obtain 

t = l - * (V) , 

therefore the general solution <j) (t,x l ) is given as an ar- 
bitrary function of the first integral X : 



(19) 



(f>(t,x*) = <S> (t + ^ (x 1 )) . 

Thus the general solution for equations (fl5|l and 1(16)1 
contains two arbitrary functions. Note that the system 
(fT7)) does not have any other general solutions besides 
(fT9")l . It is easy to prove that this solution satisfies the 
equations (fl"5")l and (fl6|) . Indeed we have 



and dj > 



d$„ _ 



therefore 



X 



~dX ) 



,oo 



[g- + 2g 0l d^ + g lk d t d k ^) 
and the time component 



(20) 



V2X y/goo + ytiQ-q, . 
along with the spatial components 



0i- 



ft* 



V2X y/goo + 2g 0l d^ + g lk d l d k ^! 

are obviously time independent because the metric is 
stationary. It is worth mentioning that by using the 
condition lfT4|) we would arrive at the general solution 
= t + * (x l ), missing the arbitrary functional 
dependence <&. Note that arbitrary field redefinitions cor- 
respond to the freedom in choosing <I>. 



C. Which Lagrangians do allow for the stationary 
configurations? 

Now let us consider the restrictions on P (<p, X) arising 
from the requirement that the pressure and energy den- 
sity should be time independent for the general solution 
((15)1 . From Eq. (TO]) we have 

d t P = + P X X = 0, (21) 

while from Eq. JTTJ 

d t (XP X ) - XP X + XP X4> <j> + XP XX X = 0. 

Eliminating X from these equations results in 

XP t x<j> - (XPxx + P,x) "5^=0. (22) 

This equation is a second order partial differential equa- 
tion for P ((f>, X). A trivial solution of this equation is a 
shift symmetric Lagrangian P (X). It is well known that 
shift symmetric theories are exactly equivalent to hydro- 
dynamics for X > 0. Obviously hydrodynamics usually 
allows for the steady flows. Let us find a general solu- 
tion of the equation ((22)l . This general solution should 
depend on two arbitrary functions. It is convenient to 
rewrite Eq. ((22)1 in the following form 



dlniPt/Px) 



d\nX 

Integrating this equation we obtain 

P,* = £r(0) XP x , 



1. 



(23) 



where a (4>) is an arbitrary function. The last equation 
([24)1 is a linear partial differential equation of the first 
order. Similarly to our previous calculations we use the 
method of characteristics to find the general solution. For 
the characteristics we have 



d<j) dX 
— — = 1 and — — 
dr dr 



(24) 



thus the integral curves are given by the equation 
The general solution of the last equation is 



X = Jexp 



a ((f)) d(f> 



where 3 is a constant of integration. Thus the general 
solution to the equations ([23)1 and f22]) is an arbitrary 
function of the first integral 3 of the dynamical system 
(O: 



P{(j>,X) = F(Xe fwS j 



(25) 



6 



where F and / {4>) = J a (<fi) d(f> are arbitrary functions. 
Note that all solutions of lf22]l are described by lf25|) . It 
is obvious that the Lagrangian ([2"5| has a hidden shift 
symmetry. Namely, we can always perform a field redef- 
inition 

4>{4>) = j d0e /w)/2 , (26) 

so that the new Lagrangian is shift symmetric P {4>, X) = 
F(X), where X = v <f> = Xe f ^\ Thus all 

scalar field theories which allow for stationary configura- 
tions are necessarily shift symmetric (explicitly or after 
field redefinition). Further, we will use the notation 4> 
always for such field variables in which the system is in- 
variant under shift transformations 4> ~ > 4> + c > where c 
is an arbitrary constant . 

Finally we can specify the profiles $ of stationary con- 
figurations. Equations lf2lj) and lf20j) yield 

P,* + Px (^) {g°° + 2g 0l d^ + g lk d t d k y) = 0, 
and using Eq. lj23j) and |20|l we obtain 

We know that in terms of the new field 4> the Lagrangian 
is shift symmetric. Thus for this parametrization <j{4>) = 
0. Therefore $ (T) — al + (3 = t + * (x i ) where we have 
absorbed the constants into ^> and t. Thus in terms of the 
field variable <j>, in which the theory is shift symmetric, 
the possible stationary configurations are always given 

by 

4> = t + * (x l ) , (28) 

and we are back to the usual ansatz (fl4| . The station- 
ary configurations in terms of the field variable 4> can be 
obtained by solving equation (J26j) or (|27|l with respect 
to <j>. This procedure determines the function <I>. While 
the function \& (a; 1 ) has to be fixed from the equations of 
motion and boundary /initial conditions. 

It is worth noting that if the metric g^ possesses an- 
other Killing vector corresponding to e.g. axial symme- 
try £eg^ v — then we can apply the result i|28p to the 
angular variable 6. Thus the solution is 

4> = t + m + * (a4) , 

where £1 is a constant and x\ denotes the rest of the 
coordinates. 

III. CONCLUSIONS AND DISCUSSION 



(may be after a field redefinition) the system has to be in- 
variant with respect to the transformation 4> — * 4> + c, for 
all constants c. The result is valid in the self-consistent 
case where the geometry is produced by the scalar field 
as well as in the test field approximation where the sta- 
tionary field configuration appears on the gravitational 
background governed by other sources. The shift sym- 
metry implies the conservation of the Noether current 

^ = P.x^^- 

Interestingly, the equation of motion implies V M J' 1 = 0, 
which is a statement of the conservation of the current 
J M . In the case when V M <^ is timelike the current J M 
can be written in the form of an effective particle density 
current J M = nu^, where the particle density 10 is 

n = V2X~P~ X . 

Note that this current is not covariant under field re- 
definitions. The conservation V M J M = of the particle 
density current usually holds in the standard hydrody- 
namics. However, the most interesting models of cos- 
mological scalar fields do not posses this additional con- 
servation law associated with the shift symmetry. Thus 
the result obtained in this paper implies that there is no 
exact Bondi (steady flow) accretion for popular classes 
of models for dynamical Dark Energy like Quintessence 
and k-Essence. This result may not have a very strong 
qualitative impact on the growth of Black Holes or on 
the evolution of the cosmological fields around them. In- 
deed, one should expect that the accretion rate should 
be in any case rather small (for the case of canonical 
scalars see Ref. [33]) • Especially in the late / current 
universe, one can almost always neglect the growth of 
the Black Hole along with the corresponding backreac- 
tion. Nevertheless, this result changes the setup for the 
investigation of the problem. Now in order to study how 
these fields could accrete onto Black Holes one is forced 
to solve the Cauchy problem for nonlinear partial differ- 
ential equations, instead of solving the boundary problem 
for nonlinear ordinary differential equations. In partic- 
ular to approach this problem one has to choose some 
initial configuration for the field and it's time derivative. 
At this stage, it is not clear what are reasonable, physi- 
cally motivated initial conditions and at what time they 
should be posed. This is very different from the case 
of Bondi accretion where the boundary conditions are 
fixed by cosmological evolution and the membrane prop- 
erty of the BH horizon. However, it may happen that 
there are some special attractor or self-similar regimes 
to which the solutions would approach in the late time 
asymptotic. Nevertheless, one cannot guarantee either 
the existence of these attractors nor their uniqueness for 



In this paper we have proved that the existence of sta- io Note that this number density is none other than the canonical 

tionary configurations requires shift symmetry. Namely momenta for the field <j> in the co-moving reference frame. 
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a general model. Moreover, even if a unique attractor 
exists, then it is not a priory known how wide the base 
of attraction is in the phase space consisting of initial 
configurations of the field and its time derivative. Thus, 
the procedure for finding these attractor solutions is not 
only a predominantly numerical exercise, but also gener- 
ically not very promising and predictive. Nevertheless, it 
is very interesting to find examples of scalar field systems 
possessing solutions of this type. In [30] it was demon- 
strated that for canonical scalars and many potentials 
the solutions indeed approach the steady flow. 

In addition one has to mention that having a shift sym- 
metric theory is a necessary, but not sufficient for the 
existence of a stationary configuration. For example, in 
hydrodynamics there can be either exceptional theories 
or even exceptional boundary conditions for which there 
are no stationary configurations. In particular the simple 
accretion of dust onto a Black Hole occurs along geodesies 
and therefore is not steady. A similar situation happens 
in the case of the Ghost Condensate for which the accre- 
tion rate blows up when the field configuration at spa- 
tial infinity ap pro aches the condensation point (compare 
Refs. |3l|and [35j]). Moreover, in the DBI model consid- 
ered in [20l. l2l| it was found that a physically meaningful 
steady state accretion is not possible when the sound 
speed at spatial infinity is c 2 s > 4/3. 

In this paper we have considered only a single self inter- 
acting scalar field. It would be interesting to study other 
types of fields, in particular one could think of scalars 



with internal degrees of freedom e.g charged scalars ac- 
creting onto a charged Black Hole. We expect that the 
appearance of new external forces and internal degrees 
of freedom can change the picture. Another interest- 
ing problem is to find possible attractor or self-similar 
asymptotic solutions and develop a perturbation theory 
around them. As we have shown stationary configura- 
tions are possible only for theories which are equivalent to 
perfect fluids. This result reveals once again that the re- 
lation between hydrodynamics and field theory is rather 
deep. Therefore we think this connection deserves a fur- 
ther study. We found that investigation of possible dy- 
namical backgrounds around Black Holes is interesting 
not only from the point of view of mathematical physics 
but may be relevant for a better understanding of both 
Black Holes physics and may be even the nature of Dark 
Energy. 
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